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ABSTRACT 

We  present  simple  inequalities  for  the  Riemann  problem  for  a  Hamilton- Jacobi 
equation  in  N  space  dimension  when  neither  the  initial  data  nor  the  Hamiltonian 
need  be  convex  (or  concave).  The  initial  data  is  globally  continuous,  affine  in  each 
orthant,  with  a  possible  jump  in  normal  derivative  across  each  coordinate  plane, 
Xi  =  0.  The  inequalities  become  equalities  wherever  a  “maxmin”  equals  a  “minmax” 
and  thus  an  exact  closed  form  solution  to  this  problem  is  then  obtained. 
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We  shall  be  concerned  with  solutions  to  the  following  differential  equation 


(H-J) 


ift  +  HiD^^)  =  0  in  72^  X  (0,  oo) 


where  HsC{R^)  and  =  (<^n , . . . ,  (^r„ )  is  the  spatial  gradient  of  if. 


We  shall  take  special  (Riemann)  initial  data.  Let  be  constants  and 

define 

u,(a;)  =  uf  if  a:,  >  0 

(1) 

u,(x)  =  Uj  if  Xi  <  0 

for  i  =  1, . . . ,  iV.  Then  take: 

N 


(2) 


foix)  =  A-\-  ^2  “t(^)  =  A  +  X  •  u(x) 


•=1 


We  wish  to  solve  H-J  with  initial  data  (2)  in  the  class  of  viscosity  solutions 
as  defined  in  [2].  The  four  properties  of  viscosity  solutions  that  we  shall  need  here 
(apart  from  existence,  uniqueness  derived  in  [2])  are: 

(Pi)  The  solution  f(x,t)  is  a  non-decreasing  function  of  the  initial  data. 

(P2)  The  partial  derivatives  f^^  satisfy  a  maximum  principle  at  points  of  continuity, 
i.e.  for  i  =  1, . . .  ,N: 


min(Ui  ,  uf )  <  <  max( ). 

(P3)  The  speed  of  propagation  is  finite. 

(P4)  If  V’(a^27  • . . ,  f)  is  a  viscosity  solution  of 


2 


for  a  constajit  ui  then 


(^(z,  t)  =  UiZi  +  ^(Z2, .  .  .  ,  z;v,  t) 
is  a  viscosity  solution  to  (H-J). 

Let 

Q  =  QxxQ2X‘--x^n 

where  for  i  =  1, . . . ,  TV: 

=  {v/  TDin(u~ ,u^)  <  V  <  niax(u“,  u^)} 


Finally  we  let,  for  i  =  1, . . . ,  iV: 

Xi  =  sign  (ut  -  uf) 

For  convenience  only,  we  order  the  indices  so  that 

Xi  —  ii  i  =  1, 2, . . . ,  _y 

Xi  = -1,  i  =  TV 

(j  might  be  0  or  iV  4- 1). 

We  now  state: 


Theorem  1.  The  viscosity  solution  .o  (H-J)  with  initial  data  (2)  satisfies: 
A -t- max  max  •••max  min  •••  min  [z  •  u  —  iff(u)] 

vsefl]  pjvefl^v 

(3)  <^{x,t)< 

A+  min  •••  min  max  •  •  •  inax[x  ■  v  —  tH{v)] 

vyeUs  Pi«Oi 


We  note  that  if  all  the  Xi  =  1,  then  this  solution  is  just  max[z  •  v  —  tH(v)],  if  all 

veil 

the  Xi  =  —1  it  is  min[z  ■v  —  tH{v)].  Otherwise  we  have  a  pointwise  inequality  which 

veil 
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gives  the  exact  solution  whenever  the  first  and  last  terms  in  (3)  are  equal.  This 
occurs  e.g.  if  H{v)  =  Hi{vi, . . .  ,Vj)  +  H2{vj+i, . . .  i.e.,  if  the  Hamiltonian 

separates,  and  in  many  other  cases. 

The  rest  of  this  paper  consists  of  the  proof  of  this  theorem,  and  some  remarks 
about  both  conservation  laws  and  numerical  approximations  to  H-J. 


It  is  easy  to  see  that  the  solution  to  the  Cauchy  problem  satisfies 

(4)  <f(x,t)  =  tg(j)-\-A  =  tg{Q  +  A 
where  g  satisfies; 

(5)  g  =  C  •  D^g  -  H{D^g)  =  -H^{gi{Q^  52(C),  •  •  • ,  9n{0) 
where  is  continuous. 

In  H-J,  we  let  r  =  t,yi  =  i,  —  Ci^  for  C  fixed,  H-J  becomes 

ifr  +  H{Dy<p)  -  (^Dyip 

(H-Ji) 

=  P>r  +  H^iDyif)  =  0 

with  the  same  initial  data  (2). 


Thus,  by  (5),  to  evaluate  g(Q  we  need  only  evaluate  —H^{Dyg)  at  y  =  0,  for 
any  t  >  0.  From  (Pi)  above  we  know  that  {Dyg)y-o  lies  in  Q  for  t  >  0.  Moreover, 
if  we  integrate  (H-J)^  from  t  =  0  to  t  =  At  we  have 

<f{(l,At)  =  A-AtH'({D,g),^a) 

(6)  =^„(0)- AtH‘(iJiVo(0).DlVo(0);CJ=¥>o(0),IllVo(0);  '  • 

Here: 


(7) 


^„(o)  = 
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where  =  {0. 0, . . . ,  1, 0, . . . ,  the  unit  vector,  and  U2  \  uJj.Uy 

is  determined  by  (6). 

This  formula  can  be  interpreted  as  a  numerical  algorithm.  Suppose  we  are 
given  a  grid 

X  j.  —  Jih^i  —  1,...,  N ,  jj  0,^1,... 

and  values  of  a  discrete  function  xf^j  =  Then  for  each  j,  we  construct  the 

piecewise  affine  function  which,  in  each  of  the  2^  orthants  centered  at  j,  interpolates 
..  xpj  and  its  N  nearest  neighbors,  for  i  =  1, . . . ,  AT.  From  (P3)  above,  if 

(8)  (CFL) 

i-l . S 

where  is  the  same  cis  0.  with  each  replaced  by  j ,  xp j ,  then  the 

solution  to  the  initial  value  problem  (H-J)*  with  the  above  affine  initial  data  in  the 
diamond  centered  at  j  when  evaluated  at  ar  =  Xj  and  ^  is  independent  of  the 
values  of  the  initial  data  outside  of  this  diamond. 

Thus  (6)  (with  <^o(0)  replaced  by  and  ip{0,At)  by  gives  us  a  mono¬ 
tone  finite  difference  scheme  approximating  which  is  in  differenced  form  with 

numerical  Hamiltonian  H^.  These  concepts  were  introduced  in  [3].  The  scheme  is 
monotone,  which  means  that  the  right  side  of  (6)  is  an  increasing  function  of  all  the 
^Pj±ei,  because  of  property  (PI).  The  function  is  called  Godrmov’s  Hamiltonian 
by  analogy  with  the  definition  of  Godunov’s  scheme  for  conservation  laws  in  one 
space  dimension  [5].  The  scheme  is  consistent,  which  means 

H^{Ui,Ui;U2,U2-,  .  .  .;Us,Un)  =  II^{Ui,U2,  . . .  ,un) 

Monotonicity  implies  that 

>  “2 ;  >  •  •  • ;  “iv) 
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is  a  nonincreEisiug  function  of  all  the  uf  and  a  non-decreasing  function  of  all  the 
u~ .  In  particular,  for  iV  =  1,  this  means  for  any  visCl  =  Qj: 

(9)  sga(u'^  -  u:;;)[H^(u^,u~)  - 

=  sgn  (u;''  -  ui)[#^(u+,u7)  -  H^{vi,Ui)] 

-t-sgn  {vi  -  uf  )  -  ^^(ui,i;x)] 

<  0 


But,  by  (P2),  H^(ut,Ui  )  =  H^(ui)  for  some  ui  in  Q.  Thus  we  have 


(10) 


^^(u+,Ui  )  =  xi  min  XiH^v^) 

vieQi 


(This  formula  was  obtained  earlier  in  [6]).  Now  we  proceed  inductively.  Suppose, 
for  iV  <  M  —  1,  we  have 


(11) 


where 


max  •  •  •  max  min  • 

••  rain  H^(vi,V2,... 

VjsQi 

,  vn) 

<  H^{uf 

^JV’  “/v) 

<  min  •  •  •  min  max 

uiefli  VjcQj  Vj+ieilj+i 

•••  max  H^(vi,V2,  . 
Vf/eSls 

■  ■ .  vn) 

Xi  =  1. 

II 

Xi  ~  1» 

II 

. 

4- 
»— » 

Next  we  have,  N  =  M  and  for  any  visQi: 

(12) 


<  0 


using  the  same  argument  as  in  (9). 
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By  the  induction  hypothesis,  this  means  we  have 

(13) 

=  X1X2  min  X2---XAf  mi?  XM  H^{vi,V2  . .  ■  vm) 

vjeCii  VMtiiu 

where  the  extrema  is  taken  on  at  V2, .  ^ .  iVm,  which  depends  on  I'l-  The  vector 
(ui,  U2, . . . ,  VN)eQi  where  uiefii  is  arbitrary.  We  next  take  min  of  the  expression 
in  (13);  If  all  the  Xi  =  1  or  Xi  =  we  have  equality  by  (P2).  Otherwise 

Xi  s  1,  I  <  i  <  j,  Xi  =  “It  j  +  1  <  i  <  M  and  we  have  the  right  hand  inequality 
in  (11).  Next  we  have,  for  any  Uj+iefij+i,  following  the  argument  above: 

(14) 

>  Vj+i ,  vj^i ; . . . ;  ,  «^) 

>  Xj+2  mm  Xj+2---XM  min  XMXi  min  Xi  •  •  •  Xj  min  x>  i’?.  •  •  • 

«j+2en,+2  VMeQm  vieJii  Vjsilj 

We  next  take  the  max  of  the  expression  in  (P4)  which  gives  us  the  left  hand  in- 

»i+i 

equality  in  (11). 


7 


I 


We  have  now  obtained  formula  (11)  for  any  N;  using  (4)  and  (5)  give  us 
Theorem  1. 

We  note  that  (11)  validates  the  conjecture  about  Godunov’s  Hamiltonian  in  [7] 
when  the  inequalities  in  (13)  and  14)  become  equalities.  That  paper  also  discusses 
the  high-order  acctirate  non-oscillatory  numerical  solution  of  (H-J)  in  some  detail. 

If  we  take  the  space  gradient  of  (H-J)  and  call  ui  =  <Pzi,U2  =  ■Pi.2,  etc.,  we 
arrive  at  the  system  of  conservation  laws 

(14)  {ui)t  +  —H(ui,...,um)  =  0,  z  =  l,...,iV 

with  intial  data; 

Ui(z,0)  =  ttj*"  if  Xi  >  0 
=  u“  if  Xj  <  0 
z  =  1, . . . ,  n 

Then  taking  the  space  gradient  of  (3)  gives  us  information  about  the  solution  to 
this  special  Riemann  problem  for  a  special  system  of  conservation  laws. 

We  finally  remark  that  if  the  initial  data  is  convex  (concave)  or  if  i3’(  ui _ _ u,v) 

is  convex  (concave)  then  the  Hopf  formulas  [1]  for  this  problem  apply  Tn  the  case 
iV  =  1  it  was  shown  in  [1]  that  these  formulas  give  the  solution  (3)  originally  derived 
in  [6].  (In  the  one  dimensional  case  the  Riemann  initial  data  is  automatically  convex 
or  concave).  The  same  must  be  true  in  the  mtJti-dimensional  convex  or  concave 
case.  Our  general  nonconvex  results  presumably  follow  from  the  rather  complicated 
formulas  in  [4],  although  the  connection  seems  to  be  unclear.  P.  Sougandis  has 
verified  for  us  (private  communication)  that  the  left  and  right  sides  of  (3)  are  always 
viscosity  sub  and  super  solutions  for  (H-J). 
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